Passive scalar diffusion as a damped wave 
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Three-dimensional turbulence simulations are used to show that the turbu- 
lent root mean square velocity is an upper bound of the speed of turbulent 
diffusion. There is a close analogy to magnetic diffusion where the maximum 
diffusion speed is the speed of light. Mathematically, this is caused by the 
inclusion of the Faraday displacement current which ensures that causality 
is obeyed. In turbulent diffusion, a term similar to the displacement current 
emerges quite naturally when the minimal tau approximation is used. Simu- 
lations confirm the presence of such a term and give a quantitative measure 
of its relative importance. 



1 Introduction 

Since the seminal paper of Prandtl (1925), turbulent diffusion has always been 
an important application of turbulence theory. By analogy with the kinetic 
theory of heat conduction, the turbulent exchange of fluid elements leads to 
an enhanced flux, T , of a passive scalar concentration that is proportional to 
the negative mean concentration gradient, 

~T = -n t VC (Fickian diffusion) , (1) 

where K t = ^u vms £ COT is a turbulent diffusion coefficient, u rms is the turbulent 
rms velocity, and £ COT is the correlation length. Equation (1) leads to a closed 
equation for the evolution of the mean concentration, C, 

This is an elliptic equation, which implies that signal propagation is instan- 
taneous and hence causality violating. For example, if the initial C profile 
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is a (5-function, it will be a gaussian at the next instant, but gaussians have 
already infinite support. 

The above formalism usually emerges when one considers the microphysics 
of the turbulent flux in the form T = u J cdt, where c w — u ■ VC is the linear 
approximation to the evolution equation for the fluctuating component of the 
concentration. Recently, Blackman & Field (2003) proposed that one should 
instead consider the expression 

dT/dt = ui + uc. (3) 

On the right hand side, the nonlinear terms in the two evolution equations for 
u and c are not omitted; they lead to triple correlations which are assumed 
to be proportional to — J^/t, where r is some relaxation time. Furthermore, 
there is a priori no reason to omit the time derivative on the left hand side 
of equation (3). It is this term which leads to the emergence of an extra time 
derivative (i.e. a 'turbulent displacement flux') in the modified 'non-Fickian' 
diffusion law, 

— — d~T 

T = -K t VC - t— (non-Fickian). (4) 
This turns the elliptic equation (2) into a damped wave equation, 

d 2 C 1 dC , 2 „ 27 , 

The maximum wave speed is obviously u rms / y/3. Note also that, after multi- 
plication with r, the coefficient on the right hand side becomes ^ru 2 ms = n t , 
and the second time derivative on the left hand side becomes unimportant in 
the limit r — > 0, or when the physical time scales are long compared with r. 



2 Validity of turbulent displacement flux and value of r 

A particularly obvious way of demonstrating the presence of the second time 
derivative is by considering a numerical experiment where C = initially. 
Equation (2) would predict that then C — at all times. But, according to 
the alternative formulation (5), this need not be true if initially dC/dt ^ 0. 
In practice, this can be achieved by arranging the initial fluctuations of c such 
that they correlate with u z . Of course, such highly correlated arrangement 
will soon disappear and hence there will be no turbulent flux in the long time 
limit. Nevertheless, at early times, (C )^ 2 (a measure of the passive scalar 
amplitude) rises from zero to a finite value; see Fig. 1. 

Closer inspection of Fig. 1 reveals that when the wavenumber of the forcing 
is sufficiently small (i.e. the size of the turbulent eddies is comparable to the 
box size), (C ) x / 2 approaches zero in an oscillatory fashion. This remarkable 
result can only be explained by the presence of the second time derivative 
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Fig. 1. Passive scalar amplitude, (C ) 1 ^ 2 , versus time (normalized by itrmsfcf) for 
two different values of ki/ki. The simulations have 256 3 meshpoints. The results are 
compared with solutions to the non-Fickian diffusion model. 




Fig. 2. Strouhal number as a function of kf/ki for different values of ReLS, i-e. 
the large scale Reynolds number. The resolution varies between 64 3 meshpoints 
(ReLS = 100) and 512 3 meshpoints (Re LS = 1000). 



term giving rise to wave-like behavior. This shows that the presence of the 
new term is actually justified. Comparison with model calculations shows that 
the non-dimensional measure of r, St = Tu rms fcf, must be around 3. (In mean- 
field theory this number is usually called Strouhal number.) This rules out the 
validity of the quasilinear (first order smoothing) approximation which would 
only be valid for St — > 0. 

Next, we consider an experiment to establish directly the value of St. We do 
this by imposing a passive scalar gradient, which leads to a steady state, and 
measuring the resulting turbulent passive scalar flux. By comparing double 
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Fig. 3. Visualizations of C on the periphery of the simulation domain at a time when 
the simulation has reached a statistically steady state, kf/ki = 1.5, ReLS = 400. 

and triple moments we can measure St quite accurately without invoking a 
fitting procedure as in the previous experiment. The result is shown in Fig. 2 
and confirms that St ~ 3 in the limit of small forcing wavenumber, fef. The 
details can be found in Brandenburg et al. (2004). A Visualization of C on the 
periphery of the simulation domain is shown in Fig. 3 for kf = 1.5. Note the 
combination of large patches (scale ~ 1/fcf) together with thin filamentary 
structures. 

Finally, we should note that equation (3) in the passive scalar problem 
was originally motivated by a corresponding expression for the electromotive 
force in dynamo theory, where the ii terms leads to the crucial nonlinearity 
of the a-effect (Blackman & Field 2002). 
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